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Abstract 

We analyze the scaling and finite-size-scaling behavior of the two-dimensional 
4-state Potts model. We find new multiplicative logarithmic corrections for the 
susceptibility, in addition to the already known ones for the specific heat. We 
also find additive logarithmic corrections to scaling, some of which are univer- 
sal. We have checked the theoretical predictions at criticality and off criticality 
by means of high-precision Monte Carlo data. 
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1 Introduction 



The two-dimensional (2D) g-state Potts model is one of the most intensively stud- 
ied systems in statistical mechanics. Although this model is very simple to formulate, 
it has a rich phase diagram. Baxter showed that for g < 4 this model undergoes 
a second-order phase transition at the self-dual point, while for g > 4 the transition 
is of first order. 

The borderline case g = 4 is the most difficult: here the transition is second-order, 
but the leading power-law scaling behavior is modified by multiplicative logarithms, 
as was first observed by Nauenberg and Scalapino and Cardy, Nauenberg and 
Scalapino These authors studied an extended Potts model (the dilute Potts 

model 0) in which, in addition to the usual g spin states, vacancies are allowed. 
The density of the vacancies is governed by a fugacity ("dilution field"). For each 
g < 4, one finds the following phase diagram: for low vacancy fugacity, there is an 
ordinary critical point belonging to the universality class of the pure Potts model; 
for high vacancy fugacity, there is a first-order transition at which the energy, the 
magnetization and the vacancy density are all discontinuous; these two transition 
curves meet at a tri critical point. In the renormalizat ion-group (RG) framework, 
there is for each g < 4 an ordinary critical fixed point, a tri critical fixed point, and a 
discontinuity fixed point (governing the first-order transition). At g = 4, the critical 
and tricritical fixed points merge, and the dilution field becomes marginal: that is the 
cause of the multiplicative logarithmic corrections. For g > 4, only the discontinuity 
fixed point survives. 

The quantitative analysis of Cardy et al. 0, |^ is based on studying the RG fiow 
for the g-state dilute Potts model, as a function of the continuous parameter g, in 
a neighborhood of g = 4. First, by using a nonlinear transformation to appropriate 
scaling fields, they obtain the RG equations in Poincare normal form through second 
order in the scaling fields and first order in e = g — 4. Next they fix the free parameters 
in this normal form by matching the RG predictions to leading order in e with the 
exactly known results for the critical exponents when g < 4 P, |^, § and for the 
latent heat when g > 4 Finally, by integrating the RG equations at g = 4 with 
suitable boundary conditions, they obtain predictions for the leading-order critical 
behavior of: a) the specific heat, the correlation length, and the magnetization as 
functions of the thermal field at zero ordering field; b) the critical magnetization as 
a function of the ordering field; and c) the critical two-point correlator as a function 
of the distance between the two points. In all these cases, the power-law dependence 
expected generically for a second-order phase transition is modified by multiplicative 
logarithmic corrections. 



^ In addition to these analytical works, the possible existence of multiplicative logarithmic correc- 
tions was proposed almost contemporaneously by Rebbi and Swendsen | ^ on the basis of numerical 
Monte Carlo renormalization group data. 

^ At that time that Refs. H, ^ were written, the exact formulae for the critical exponents were 
only conjectured. They have now been confirmed by both Coulomb-gas and conformal-field-theory 
methods: see Appendix A.l for a brief review. 



2 



In this paper we show that by extending the analysis of Cardy et al. through 
third order in the fields, we can obtain not only the leading critical behavior but 
also the universal leading corrections to scaling: these are additive terms of the 
generic form log log / log, with universal amplitudes that we can compute exactly. In 
addition, there are nonuniversal subleading correction-to-scaling terms of order 1/ log. 
A similar behavior arises in other models with marginally irrelevant operators, such 
as the four- dimensional cp'^ (or Ising) model |T^, the three-dimensional tricritical 
ip^ (or spin-1 Ising) model |jTl|, and the three-dimensional ferromagnet with strong 



dipolar interactions |T^. 

A second goal of this paper is to extend the renormalization-group analysis of 
Cardy et al. to obtain the finite-size-scaling (FSS) behavior, as well as the leading 
corrections to finite-size scaling, for the 4-state Potts model in a periodic L x L box. 
Finite-size scaling [|T4| , pTS] , |16|] plays an important role in the analysis of Monte Carlo 
simulations, which of course deal with finite systems and must thus be extrapolated to 



the infinite-volume limit. This extrapolation [0, |T8|, |T^] is a very delicate procedure, 
and the corrections to finite-size scaling induce systematic errors in the extrapolation. 

For q < 4 the FSS behavior of the Potts model is well understood. The critical 
singularities are rounded-off: at finite volume, the specific heat (resp. the suscepti- 
bility) has a peak of height ~ L"/'' (resp. Li^/'^)^ and this peak is shifted from the 
critical point by an amount of order L~^^'^. The corrections to this leading behavior 
are suppressed by a factor L~'^, where u is a correction-to-scaling exponent. 

For g > 4 the FSS behavior is that appropriate to a first-order phase transition: 
the height of the peak in the specific heat scales as ~ L'^, and the shift as ~ L~'^, 
where d is the dimensionality of the lattice. The corrections are of order and L"^"^, 
respectively. This behavior has been proven rigorously for g ^ 1 by Borgs, Kotecky 



and Miracle- Sole [pO] , |21| , P2| , and has been confirmed numerically by Billoire et al. 
^ for g = 20. 

For q = 4 the FSS behavior is more complicated, but it can be obtained by a 
relatively trivial extension of the infinite-volume RG analysis. Not surprisingly, we 
find multiplicative logarithmic corrections to the leading power-law behavior (some 
of which have been found previously p4|); in addition, we find additive corrections of 
orders log log L / log L and 1 / log L. In all cases we try to distinguish which corrections 
are universal and which are nonuniversal. 

We remark that the FSS behavior of the eigenvalues of the transfer matrix in 



an L X oo strip is also known [^, ^: there are additive logarithmic corrections to 
scaling due to the marginal dilution field. 

The third goal of this paper is to test the predicted finite-size-scaling behavior 
for the 4-state Potts model by means of a high-precision Monte Carlo simulation 
using a Swendsen- Wang-type algorithm. We find that our Monte Carlo data are 
consistent with theory if we include both the predicted multiplicative logarithms and 
the additive logarithmic corrections. However, it would have been impossible to 
deduce the theory from these data, which are also consistent with other functional 
forms. In particular, we are unable to observe numerically in an unequivocal way the 
correct powers on the multiplicative logarithms or the correct universal amplitudes on 
the additive log log L/ log L corrections: these features are obscured by the presence 
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of the nonuniversal 1 / log L corrections. 

The plan of this paper is as follows: In Section |^ we carry out the RG analysis for 
the 4-state Potts model in infinite volume, with emphasis on the multiplicative loga- 
rithms and on the universal additive logarithmic corrections. In Section ^ we extend 
this analysis to determine the finite-size-scaling behavior and the leading corrections 
to it. In Section ^ we describe our Monte Carlo simulations. In Section ^ we compare 
our numerical data, both on and off criticality, to the RG predictions. In Section ^ we 
report briefly our data on the dynamic critical behavior of our Swendsen- Wang-type 
algorithm. In the Appendix we show how the Cardy et al. RG analysis for q A can 
be extended to cubic order in the fields, yielding the universal third-order term in the 
RG flow for the dilution field. 



2 Scaling equations in infinite volume 
2.1 Renormalization-group flow 



It is well known [g, y, ^ that the scaling behavior of the two-dimensional 4-state 
Potts model has logarithmic corrections. This is due to the presence in this model of 
a marginal operator, which is absent in any other two-dimensional Potts model. 

Cardy, Nauenberg and Scalapino [|, H described the renormalization-group flow 
for the g-state dilute Potts model by using three scaling fields: a thermal field an 
ordering field and a dilution field ip^ The pure Potts model corresponds to a large 
negative value of the dilution field. Near criticality, the thermal field is proportional 
to the temperature deviation from criticality (0 ~ Jc — J). The fields and h are 
relevant, and the field is marginal when g = 4. The critical point is found at 
= /i = 0. In the larger space {q, (j), h,ip), the point q = 4:, (j) = h = ip = 
is a multicritical point, at which a curve of critical points [q < 4, (p = h = 0, 
ip ~ — (4 — g)^''^] meets a curve of tricritical points [g < 4, = /i = 0, ■i/' ~ +(4 — g)^/^]. 

Cardy et al. found that the renormalization-group equations for g = 4 under an 
infinitesimal change of scale dl, keeping terms through second order in the fields, are 



#(/) 

dl 
d(t){l) 

dl 
dhjl) 

dl 

where 



ai){lf (2.1a) 
[yr + brnW) (2.1b) 
[yH + cmni) (2.1c) 



3 15 1,3 1 , , 

yT = -, yH = —, a = -, b = —, c=— . 2.2 

2 8 TT 47r lOTT 



^ Of course, in addition to these three fields there will be an infinite number of irrelevant scaling 
fields. But these will produce corrections to scaling that are suppressed by powers of the length scale, 
hence are negligible compared to the logarithmic corrections arising from the marginal operator tp. 
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The solution of these equations is: 



^(n = "^^^^ (2.3a) 



where I is the logarithm of the length-rescaling factor. From (|2.3 b,c) it is clear 
that the thermal and ordering fields are relevant (i.e., they grow exponentially as 
/ oo). The dilution field ip is marginally irrelevant when '?/'(0) < 0, in the sense 
that ipil) — > as / — >■ cxD but at a subexponential rate. This induces multiplicative 
logarithmic corrections (i.e., powers of I) in and h\ only when ip{Q) = do we have a 
pure power-law (i.e., exponential-of-/) behavior for the two relevant fields. This latter 



behavior occurs in the Baxter- Wu model [271, which is believed to be in the same 



universality class as the 4-state Potts model but does not exhibit any multiplicative 
logarithmic correction. We shall henceforth assume without comment that < 0. 

The equations ( |2.1| ) are the Poincare normal form of the RG flow at g = 4, 
taking into account terms through second order in the fields 0]. The analysis of these 
equations yields the correct leading terms (power and multiplicative logarithm) for all 
the critical observables (correlation length, specific heat, etc.). However, we should 
include more terms in order to be able to analyze the corrections to scaling to these 
quantities. Equation (p.l|a) should be replaced by 

^ = F{m) = am' + a'm' + ■■■ m 

where a' is some constant, and the dots • • • stand for higher powers of Note that 
the ratio a'/a^ is universal, in the sense that it cannot be altered by any smooth change 
of variable; the coefficients of order and higher, by contrast, can be set to any 
desired values (e.g. zero) by a smooth change of variable ip ^ ip + a2'ip'^ + a^ip'^ + ■ ■ ■ . 
The solution of (|2.4|) has the following asymptotic behavior for large /: 



^(0 = 

al 



(2.5) 



Note that the log/// correction term here is universal: it does not depend on ip{Q) 
or on the parametrization of ip- By contrast, the 1// corrections are nonuniversal, as 
their value depends explicitly on V'(O). The terms of order and higher in (|2.4|) , 
even if present, affect only the 1// corrections. 

The effect of the term a'ip{l)^ in ( p. 41) on the solutions of (|2.1| b,c) can be computed 
easily. Let us define the function /a,b(/; '?A(0)) to be the ratio 9{l)/9{0) when solving 
the generic equation 

^ = [A + Bam] 0{l) , (2.6) 
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where ■?/'(/) is the exact solution of ( p.4|) with the specified initial value Namely, 
the function fA,B is given by the integral 



exp (^JjA + Baij{l')]dl'^ 



exp 



^(0 

\A + Bail)', ^, , , 
^(0)^ ^'F{i,') 



(2.7a) 
(2.7b) 



The asymptotic behavior of j a,b for large / is 



hAhm)) = e-^' (-ai/.(0)0" 



(2i 



Thus, the solutions (p73|b,c) are replaced by 

0(O)e3'/2 



HI) = h{0)f^.{l;m) 



(-aV^(0)/)3/4 
/i(0)e^5V8 
' (-aV'(0)/)Vi6 



3/4 



(2.9a) 



-1 1/16 



(2.9b) 



One must also include higher powers of ■?/'(/) inside the square brackets in (p.l|b,c), but 
these terms cause only subleading corrections, multiplying (2.9) by factors 1 + 0(1//). 

The value of a' can be obtained by studying the RG fiow in a neighborhood of the 
multicritical point g = gc = 4, and matching the exponents with the exactly known 
value of the next-to-leading thermal exponent as a function of q. The details of 
this computation are given in Appendix 0; the result is 



27r2 ■ 



(2.10) 



2.2 Correlation length 



In infinite volume one can define various different correlation lengths: the most 
important of these are the exponential correlation length (= inverse mass gap) 



(exp) 
oo 



\X\ 



lim - 

|a;|^oo logG(Xj 



(2.11) 



(where x is taken to infinity along a coordinate axis) and the second-moment corre- 
lation length 



. E x ^G{x) 

cm = 1 1_ jn 

' 2d EG(x) 



1/2 



(2.12) 



where G{x) is the two-point correlation function and d is the spatial dimension. 
(We have appended a subscript oo to emphasize that these quantities are defined in 
infinite volume, in order to distinguish them from the finite-volume quantities to be 
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introduced in the next section.) Both of these correlation lengths, as well as other 
similar ones, are expected to scale in the same way near the critical point (but with 
different prefactors). 

Let, therefore, ^oo be any one of these correlation lengths. It behaves under a 
change of scale / as 

u{m,hio),m) = e'u{m.Ki)M))- (2.13) 

[Henceforth we shall drop the arguments on the initial values of the fields, and 
simply write (p = 0(0), h = h{0), ip = ip{0)-] Substituting the solutions (2.9), we get 
the scaling form 

U{<P,h,i;) = e'uUhMl;i;), i (/;V^), m) (2-14) 



where is given by 

Let us now define I* = l*{(f>, ip) as the solution of the equation 

(This solution is unique when |0| is sufficiently small.) The exact form of /* as a 
function of and ip cannot be obtained in closed form. However, we can get the 
following asymptotic expansion for small (p (and fixed ip < 0): 

/*(0,^) = log 101 + ^log(-log|0|) + llogf^] 



(3 3a' \ log(-log|^,, . 

Note that the first, second and fourth terms in (|2.16|) are universal, while the third 
term and the 0(1/ log |0|) correction are nonuniversal (they depend on ip). It follows 
from (|2.16|) that for any exponents A, B we can write 

F^,b(0;V^) = /^,B(r(0,0');V') (2.17a) 



-2^/3^- log i0n(^-^) 



A o\ I —2aip\ 



X 



3 ; 



4 2a2 / V2 / -log 101 \log\ 



(2.17b) 



Note that the leading power and log are universal, as is the additive correction pro- 
portional to log(— log 101)/ log 101; the constant factor and the 1/ log |0| additive cor- 
rection are nonuniversal. 

Let us now insert / = l*{(p,ip) into ( |2.14| ): using ( ^.l(j| ) / (2.17), and restricting for 
simplicity to zero field {h = 0), the result is 

eoo(0,O,^) = Fi,o(0;^) eoof±l, 0, 2alog|0| ^^'^^''^ 
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1-2/3 



(2aiij 



1/2 



log 



V 3 

xeoo(±i, 



Is 4a2i -log I 



+ 01 



log 



2a log I 



+ 



(2.18b) 



1-2/3 



log 



1/2 



'3 _ 3a^\ log(-log|(/>|) 
8 4a2] -log 101 



log 



where 



1/2 



eoo(±l,0,0) 



(2.18c) 



(2.19) 



is a nonuniversal amplitude; here the + (resp. — ) sign corresponds to the high- 
temperature (resp. low-temperature) side of criticality. In the passage from ( p.l8| b) to 
( ^.18| c) we have assumed that ^oo(±l, 0, 0) 7^ and that ^oo(±l5 0, ip) is a smooth func- 
tion of ip near = 0: together these assumptions imply that ^oo(±l, 0, 3/(2a log |0|) + 
• • ■) = const + 0(1/ log 101). Physically, this amounts to assuming that -0 is a non- 
dangerous irrelevant variable |3y, |31]]. We shall henceforth make this assumption 
without further comment. The leading power and multiplicative logarithm in ( 2.18| c) 
were first obtained in the universal log(— log |0|)/ log |0| additive correction is 
new. 

If we invert ( ^.l^ c) we get 



1 - 



3a' log log 
4a2 log ^00 



+ 01 



log^o 



where 



\3/2 



(2.20) 



(2.21) 



is a nonuniversal amplitude. Again, the correction of order log log ^00/ log ^00 is uni- 
versal, while the correction of order 1/ log ^00 is nonuniversal. 

To study the critical isotherm (with h > 0) we first define the scale /** = l**{h, ip) 
to be the solution of the equation 



This solution is given asymptotically for small h by 



r{h,ij) = -— log/i + — log(-log/i) 
15 3U 



(2.22) 



— log 

30 ^ 



-2a^\ 



480 



16a2 



log (-log h) 
— \ogh 



+ O 



15 ; 

(— 



+ 



(2.23) 
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We can also define a function analogous to Fa,b- 

= h-'^/''{-logh)^o-B 



[2.24a) 



-B 



15 



X 



lie 8a2 J V30 / -logh Uogh 



(2.24b) 



If we now choose / = l**{h,i)) in (|1D, we get on the critical isotherm (0 = 0) 



eoo(0,/i,V^) = Gi,o(/i;V^)eoo(^0, 1, 

^ /,-8A5(_log/,)l/30 



15 1 



+ 



8a log h 

( 1 a' \ log(— log/i) 



(2.25a) 



480 16a2 



log /i 



O 



1 



log /i^ 
(2.25b) 



2.3 The free energy and its derivatives 

The singular piece of the free energy per unit volume behaves under a change of 
scale / as 

/sing(0(O),/i(O),^(O)) = e-2%„g(0(Z),/i(/),^(O). (2.26) 

If we plug the solution (2.9) of the renormalization-group equations (|2.1| ) into ( |2.26| ), 
we get 

/,ing(0,/i,V^) = e-273ing(0/|,|(/;^), /^/f ,^(^;^), m) (2.27) 

where ip{l) is given by ( p.5|) . 

Let us now choose the scale / = l^i^cpjip) as in ( |2.15| ); using ( p.l6| )/(2.17), we get 



2a log 



+ 



(2.28) 



where the + (resp. — ) sign corresponds to the high-temperature (resp. low-temperature) 
side of criticality. At zero field {h = 0) this expression can be written as 



/3ing(0,O,^) ^ D± 101^/3 (-log 
where 



1 _ [ 3 _ 3a^\ log(-log|0|) 

1 4 2a2] -log 101 llogi: 



-2a^^ 



/sing(±l,0,0) 



(2.29) 



(2.30) 



is a nonuniversal amplitude. The leading power and multiplicative logarithm in ( p. 291) 
were first obtained in 0]; the universal log(— log |0|)/ log \(f)\ additive correction is new. 
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If we differentiate (|2.28|) twice with respect to 0, we get the specific heat. At zero 
field {h = 0) we get 



Ch(0,O,^) 



1 _ _ 3a^\ log(-log 
1 4 2a- -logic/ 



01 



log 



(2.31) 



The leading behavior of the specific heat was previously obtained in 

If we differentiate ( |2.28| ) once with respect to the ordering field h, we get the 
magnetization. The result at zero ordering field {h | 0) in the low-temperature 
regime (0 < 0) is 



M(</.,0,^) 



8 ' 16 8 ' 16 

(_0)l/12(_log(_0))-l/8 X 

'3 3a' \ log(-log(-0|)) , 



1 



32 16a- I -log(-^ 



01 



log(-^ 



(2.32a) 



.(2.32b) 



If we differentiate ( 2.28 ) twice with respect to the ordering field h, we get the suscep- 
tibility. The result at /i = is 



X(0,O,^) 



F_2,o(0;^)F^ !(</>; ^)' = F7 ^U;^p) 

8 ' 16 4 ' 8 

\<P\-'/%-\og\<P\f/' X 

'9 9a' \ log(- log 
16 ~ 8a2 J -log|^ 



1 - 



+ 



log 



f2.33a) 



. (2.33b) 



The leading power and multiplicative logarithms in (2.32b) were first obtained in [Q. 
The result (2.33b) did not appear in 0, but it can of course be obtained directly 
from their approach. The universal log(— log |0|)/ log |(/)| additive corrections are all 
new. 

Finally, we can write the asymptotic behavior of the specific heat, the magnetiza- 
tion, and the susceptibility as functions of ^oo- Using ( |2.2CI| ) we get 



C 



H 



M 



X 



-1/8 



f7/4 
(logeoo)^/^ 



1 

1 + 
1 + 



3a' log log 
2a2 log ^oo 

a' log log 
16a2 log ^oo 

a' log log 
8a2 log ^oo 



O 



1 



log^o 



O 



+ 01 



Jog^oo 

1 \' 



log^o 



(2.34a) 
(2.34b) 
(2.34c) 



Note that the (universal) corrections of order log log^oo/ log^oo are proportional to 
the ratio a' /a- and to the exponent of the multiplicative logarithm, since the other 
terms of the same order cancel out. 
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Finally, we observe that the following two relations hold: 



OO T 



X 



1 + 
1 + 



1 



.log I 
1 



log 



f2.35a) 



(2.35b) 



Note that in these relations there are neither multiplicative logarithmic corrections 
nor additive log(— log |0|)/log|0| corrections. It follows that the hyperscaling laws 



du 



2- a 

7 + 2/3 



(2.36a) 
(2.36b) 



are satisfied without logarithmic corrections. 



The magnetization on the critical isotherm (0 = 0, > 0) can be calculated by 
first choosing the scale / = l**{h,ip) as in (g]^. Plugging (P^)-(2.24) into { ^^ ) 
we get 



/smg(0, h, t/j) = G^2,o{h; tp) Ung 0, 1, 



15 



V"' " 8a log ^ 
/ii6A5(_iog/,)-i/i5 X 



+ ■ 



^ f 1 a' Yogi-\ogh) ^ f 1 



(2.37a) 



. (2.37b) 



^240 So?) -\ogh ' \\ogh^ 
By differentiating this equation with respect to h we get the critical magnetization: 

M(0, h, ^) ~ h^/^\- log h)~^^^^ 



1 o' \ log(- \ogh) ^^f 1 



240 8a2 



log h 



\ogh 



(2.38) 

The leading power and multiplicative logarithm in this result were first obtained also 
in 0. The universal log(— log /;,)/ log /i correction is new. 

2.4 Two-point correlators 

The authors of also included the following behavior for the two-point correlation 
function at criticality: 



G{x; 0,0, ip) 



1 



|x|V4(log|x|)V8 

If we assume a natural scaling law of the form 

1 



as X — * OO . 



(2.39) 



G{x; 0,0, z^) 



\x\ 



|a;|i/4(log 1x1)1/8 ^«U^(0, 0,7/;) 



(2.40) 
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where Fq is some scaling function, then we can deduce equation ( |2.34c|) through the 
integral 

X = [(fxG{x). (2.41) 



Moreover, we can deduce the relation (2.34b) by taking |x| ~ .^oo and using the 
hyperscaling relation 

Gd^l-^oo) ~ lim G{x) = M\ (2.42) 

|a;|^oo 

3 Finite-size scaling 

In this section we analyze the finite-size-scaling behavior of the two-dimensional 
4-state Potts model on an L x L lattice with periodic boundary conditions. 

3.1 Renormalization-group flow in finite volume 

To obtain the finite-size scaling of this model we only have to adjoin a new "scaling 
field" L~^, where L is the linear size of the system. The fixed-point value of this field 
obviously is = 0, corresponding to the infinite- volume limit. The behavior of this 
scaling field under a change of scale is trivial: to (P?T|a-c) we need to adjoin the fiow 

''^"'W=L-(0. (3.1) 



dl 



3.2 Correlation length 



Not all of the correlation lengths employed in infinite volume have sensible ana- 
logues in a fully finite lattice: for example, the exponential correlation length ( ^.111 ) 
makes sense only if the lattice is infinite in at least one direction. However, the 
second- moment correlation length ( 2.12 ) does have a sensible extension to finite vol- 



ume; this extension is, however, not unique. One reasonable definition for a periodic 
lattice of linear size L is 

^ - 2sin(7r/L) ' ^^'^^ 

where x is the susceptibility (i.e., the Fourier-transformed two-point correlation func- 
tion at zero momentum) and F is the corresponding quantity at the smallest nonzero 
momentum (27r/L,0): 

X = Y.G{x-L) (3.3) 

X 

F = ^e2^'^^i/^G(x;L) (3.4) 



X 
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where G{x; L) is, of course, the two-point correlation function on the finite lattice. 
Another definition could be 



V 



1/2 












27r 





1/2 



(3.5) 



Hereafter, ^ will denote any reasonable finite-volume correlation length. 
The generalization of (|2]g)/(|2]T|) is 

a<l>,h,^,L-') = e'mi).h{l)Ml),L-\l)) 

y 2'4 8 '16 J 

where ■?/'(/) is given by ( p.5|) . Let us now choose the scale 

/ = log L , 

so that (3.6) becomes 

a<P,h,i;,L-') = Le(0/|,3(logL;V^), /i/^,^(logL;V^), V^(logL), 1 



(3.6a) 
(3.6b) 

(3.7) 
(3.8) 



One quantity of interest is the L-dependence of the correlation length at criticality 
= h = 0). From ( |3.8D and ( ^.51 ) it is trivial to see that 



e(0,0,V^,L-i) = L^(0, 0, V'(logL),l) 



. A ^ log log L ^1 
X + r + B + O 



logL 



{log ly 



where 



X* = e(0, 0,0,1) 
1 



A 
B 



a dz 



A 



(0,0, z,i; 



z=0 



(3.9a) 
(3.9b) 

(3.10a) 
(3.10b) 

(3.10c) 



are all universal; the corrections of order l/(logL)^ are nonuniversal. 

We can also obtain a finite-size-scaling law off criticality. For simplicity let us 
restrict attention to /i = 0. From (13.81) we know that at zero field 



d(/./i3(logL;V'), 0, 



1 



alogL 



+ 



(3.11) 



And from ( p.8|) we know that 

0/3 3 (log L; 7/-) = 0L3/2(-a^logL)-3/^ 



i + ^i^ + ol 

log L ' 



logL^ 



(3.12) 
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where the 0(1/ log L) correction term is nonuniversaL Therefore, we obtain the FSS 
law 



^((/.,0,V;,L-i) 
L 



-3/4N 



+ o 



log log L 



(3.13) 



where the function is universal modulo a nonuniversal rescaling of its argument. 

The disadvantage of this approach is that the corrections to FSS are of order 
log log L/ log L. A better way to write the FSS law is to consider the relation between 
the finite-volume and infinite- volume correlation lengths at the same 0, h, ip. Again 
restricting for simplicity to h = 0, let us insert in (|3.11|) the expression (|2.20|) for (p 
as a function of ,^oo = ^oo(05 0, ip). The first argument on the right-hand side of ( ^.111 ) 
can be written (for L,C,oo ^ 1) as 



(/./,,(logL;^) = ±E± 



-3/2 /, ^ 



L 



logL 



I ^ ^ /loglogL _ log log^o 
4a2 I logL log^oo 



3/4 



A' 

+ r + 



B' 



logL log^o 



+ 



(3.14) 



where E± = ^^,(±1,0,0)3/2 

is a universal amplitude, the coefficients A' and B' are 
nonuniversal (they depend on tp) , and the dots denote corrections of order log log L/ (log L)^ 
and loglog^oo/(log^oo)^- If "we denote hj r] = ri{ip) the coefficient of the nonuniversal 
l/l correction in the formula ( p.8|) for then the values of A' and B' are given by 



(3.15a) 
(3.15b) 



A' = -^rji^) 

B' = -^r/(V^)-^logeoo(±l,0,0)-^^logeoo(±l,0,x) 



x=0 



If we now consider the finite-size-scaling limit 
fixed, we get 



0, ip fixed, L — s> oo with ^oo/L 



>/3 3(logL;V^) = ±E± 



-3/2 



1 + 



A' + B' + \\og{U/L) 
logL 



+ 



(3.16) 



Miraculously, the combination A' + B' appearing here does not depend on ip. There- 
fore, inserting this into ( p.ll[ ), we obtain 



'eoo(0,o,^)^ 



+ F, 



'eoo(0,O,^)^ 



logL 



+ 



(3.17) 



where the functions and F^ are universal (they do not contain any piece depending 
on ip). It is quite remarkable that not only the leading finite-size-scaling function, but 
also the leading correction to it, is universal (in amplitude as well as in shape). This 
is a special feature of the logarithmic corrections to scaling induced by a marginally 
irrelevant operator, and is not observed in the more common context of power-law 
corrections to scaling induced by irrelevant operators. 
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The equation ( |3.17|) can be formally inverted, yielding 



L 



Again both functions and are universal. 



3.3 The free energy and its derivatives 

The finite-volume generalization of ( |2.271 ) is 



(3.19) 



where ■?/'(/) is given by ( p.5|) . As in (|3^ ) we choose / = logL, and get 

/sing(0,/i,^,^"^) = L"Vsing ^0/3,1 (log kf is ^^{\og L] Ip) , - ^ ^ + " ■■, 1^ ■ 

(3.20) 

The finite-size behavior near criticality of the specific heat, magnetization, and sus- 
ceptibility can be obtained by performing the appropriate derivatives with respect to 
(j) or h. For simplicity let us restrict attention to zero field, where we ge 



C^h(0,O,7/;,L 



(-aV^logL)3/2 



^ ^ 3a' log log L ^ kcjj ^ 
2a2 log L log L 



J sing 



X, 0, 



1 



alogL 



+ ■ 



a;=</>/3 3 (log-J-;V') 



(3.21) 



x{4>,o,'4j,l 



^— aV'logL)V8 



1 + 



a' log log L k 



Sa^ log L 



+ 



sing 



(^/a 3(logL;^), y. 



logL 
1 



alogL 



+ ■ 



+ 



(3.22) 



We refrain from writing also 



M(0,O,V,i-') = 



(-ai/; log L) 1/16 

^/sing 



(9y 



1 _|_ log log £ _^ fcA/ _^ 
16a^ logL logi 
1 



(/-/i i(logL; V), y 



alogL 



l;=0 



because the finite- volume magnetization of course vanishes at /i = 0, for all (j) (by virtue of the Potts 
symmetry). This scaling is indeed valid, but the prefactor {dfsing/dh){4>, h, ■0, 1) vanishes at /i = 0. 

The scaling L~^/^{[ogL)~^/^^ at criticality applies not to the usual magnetization M = L~'^{A4), 
but rather to the "absolute magnetization" M = L^'^{\A4\). These two quantities are essentially 
identical in the low-temperature phase (/i | at fixed < 0), but are very different in the critical 
regime. Unfortunately, much of the literature (especially numerical work) has sloughed over the 
distinction between M and M. 
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For future convenience, we have denoted by ko the coefficient of the nonuniversal 
1/logL correction term within the square brackets, which arises from the nonuni- 
versal l/l term in ( p.SI ). In addition to the two correction terms within the square 
brackets, we have (through order 1 / log L) three other corrections arising from the 
factors involving /sing: a universal 1/logL term arising from the derivative of the 
relevant scaling function with respect to ip ai ip = 0; and a universal log log L/ log L 
term plus a nonuniversal 1/logL term, both arising from the fact that the ffist ar- 
gument of /sing is not exactly proportional to 0L^/^(log L)~^/^, but is rather given by 
( p.l2| ). The simplest case is at criticality (0 = 0): then the second and third terms 
vanish, and the ffist one gets amalgamated with the nonuniversal terms ko/ log L. 
We thus get 



(log L) 3/2 
(logL)V8 



^ ^ 3a' log log L ^ kci^ ^ _ _ 
2a2 log L log L 

a' log log L k' 

1 H — \ ^ + ■ ■ 

Sa^ log L log L 



(3.23) 
(3.24) 



where k'r< and k' are nonuniversal. 

In conclusion, all the observables which are computed through derivatives of the 
free energy have at criticality a multiplicative logarithmic piece, proportional to some 
power of logL. They also have additive corrections of order log log L/ log L, the 
coefficient of which is —a'/a^ times the exponent of the multiplicative logarithm. 
These corrections arise together from the fact that with each derivative we gain a 
factor /^^ (log L; '?/'). On the other hand, the 1/logL corrections in ( ^.23] )/ ( ^.24] ) 
are nonuniversal because the corresponding terms in /^ ^(logL; ■?/') are nonuniversal. 
Both the multiplicative logarithm and the additive log log L/ log L term are absent 
in the correlation length (3.9). The leading power and multiplicative logarithm in 
( p.23| ) were obtained previously in [|^. It is noteworthy to remark that in some of 
the literature it was (wrongly!) assumed that the leading term for the susceptibility 
had no such logarithmic corrections [0, 0, |3^ . 

Off criticality, we get FSS laws 



X(0,O,^,L 



(logL)3/2 

(log L) 1/8 



F,,(0L3/^(^logL)-3/^)+o(W^ 
F,(0L3/^(^logL)-3/^) + o(i^)" 



(3.25) 



(3.26) 



The disadvantage of ( p.25| ) / ( p.26D is that the corrections to FSS are of order log log L/ log L. 
A better way to write the FSS law is to use C,oo/L as the argument of the scaling func- 
tion. Using the relation ( |3.16| ), we get 



Ch(0,O,7/',L- 



L 3a' log log L kc„ 

(logL)3/2 [ +2^ logL ^k^^" 
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LV4 



(log L) 1/8 



1 + 



L J "\ L 

a' log log: L 



logL 



+ 



+ 



8a2 log L log L 



+ ■ ■ 



'eoo(0,o,v^)^ 



X 

1 



logL 



+ ■ 



(3.27) 



(3.28) 



where the scaling functions Fq and Fq are universal and can be expressed in terms of 
derivatives of /sing- Thus, the nonuniversal terms kol^ogL appear only in the prefac- 
tor in square brackets. There is, of course, also a nonuniversal prefactor multiplying 
everything. 

Finally, we can re-express ^oo/L as a function of ^(0, L)/L = ^{(f),0,ilj, L~^)/L, 
using (|3.18|) , and get 



L 



(logL)3/2 



^ ^ 3a' log log L ^ kcjj ^ _ 



2a2 logL 



logL 



logL 



(3.29) 



(logL)V 



a' log log L ky 

1 H — ;t r ; — h :r-^ + ■ 



^a^ log L 



logL 



logL 



+ 



(3.30) 



The scaling functions Go and Gq are universal. 

If we study the ratio CH{4>y^y'4'y L~^)/Ch{0,0,iP, L~^), we can see easily from 
( p.29| ) that the prefactor in square brackets cancels out, and we have 



G^(0,O,V;,L"^) 
Gj^(0,0,^,L-i) 



^Ch T + ^C, 



(1) (a^.LY 



logL 



+ 



(3.31) 



where both scaling functions are universal. Therefore, both the universal log log L/ log L 
corrections and the nonuniversal 1 / log L corrections to the specific heat Ch{<Pi 0, L~^) 
come from its value at criticality. Similar reasoning applies to the susceptibility. 



4 Description of the simulations 



4.1 The Monte Carlo algorithm 

The Monte Carlo (MC) algorithm we used in our simulations was actually an 



algorithm ||3^, ^ to simulate the Ashkin- Teller (AT) model ^ . This model is a 
generalization of the Ising model to a four-state model, and it includes as a particular 
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case the 4-state Potts model. The general AT model assigns to each lattice site x two 
Ising spins = ^1 and = ±1, and they interact through the Hamiltonian 

-f^AT = -J^cr^ay - J' ^T^Ty - K^a^T^ayTy , (4.1) 

{xy) (xy) (xy) 

where the sums run over nearest-neighbor pairs (xy). The line J = J' = K is the 
4-state Potts model with /3 = 4J: 

-f^Potts = -J'^icr^ay + T^Ty + a^j^r^cTyTy) = -4:J^Sa,^aySr^,ry + const . (4.2) 

{xy) (xy) 

The critical point of the square-lattice 4-state Potts model lies at J = J' = = 
J log 3. The plane = of the general AT model corresponds to two non-interacting 
Ising models with nearest-neighbor constants J and J', respectively. 

Wiseman and Domany |^ proposed an algorithm of Swendsen-Wang (SW) type 



for the general AT model. This algorithm (called the "direct algorithm" in |3^, |34 
reduces to the standard SW algorithm ||3^ at the Ising and 4-state Potts subspaces. 



In |3J] we proposed an embedding variant of their algorithm; it reduces to the 
standard SW algorithm at the Ising plane, but not at the 4-state Potts line. How- 
ever, we gave numerical evidence that the two algorithms lie in the same dynamic 
universality class. In particular, we found that at criticality 



direct 



mt,£ 



—embedd 
'int,^ 



1.516 ±0.035, (4.3) 



where Tint.f denotes the integrated autocorrelation time of the energy (this is roughly 
the slowest mode in SW-type dynamics). However, because the embedding algorithm 
requires 1.9 times as much CPU time per iteration as the direct algorithm (since in the 
AT formulation there are twice as many spin variables as in the Potts formulation), 
our algorithm is about 25% less efficient than the standard SW algorithm. 

Let us review briefly our embedding algorithm (more details can be found in ||34|| ) . 
First, consider the Boltzmann weight of a given bond (xy), conditional on the {r} 
configuration (i.e., the r spins are kept fixed): it is 

W^bond(a., ay- r,, Ty) = e~^'^^^^-^y^ + [l - e-^'^^^^^^y^] 5.,,.^ . (4.4) 

We can simulate this system of a spins using a standard SW algorithm. The effective 
nearest-neighbor coupling 

Jfy = J{l + T,Ty) (4.5) 

is no longer translation- invariant, but this does not matter. The key point is that the 
effective coupling is always ferromagnetic. An exactly analogous argument applies to 
the {r} spins when the {cr} spins are held fixed. 

The embedding algorithm for the 4-state Potts model has therefore two parts: 

Step 1: Update of {cr} spins. Given the {r} configuration (which we hold 
fixed), we perform a standard SW iteration on the a spins. The probability pxy 
arising in the SW algorithm takes the value p^y = 1 — exp[— 2J(1 -|- T^Ty)]. 
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Step 2: Update of {r} spins. Given tlie {a} configuration (wliicli we liold 
fixed), we perform a standard SW iteration on tlie r spins. Tlie probability p^y 
arising in the SW algorithm takes the value p^y = 1 — exp[— 2J(1 + (JxCTy)]. 

One iteration of the embedding algorithm consists, by definition, of a single ap- 
pheation of Step 1 followed by a single apphcation of Step 2. 



4.2 Observables to be measured 



Let us begin by defining some basic observables. The observables of interest 
involving only the a spins are 



Ea — ^ ] ^x^y 



1 



2i7Txi/L 



+ 



2inX2/L 



(4.6) 
(4.7) 

(4.8) 



where L is the linear size of the system (we always use periodic boundary conditions) 
and {xi,X2) are the Cartesian coordinates of the point x. The observable JF^^ can be 
also seen as the square of the Fourier transform of a at the smallest allowed non-zero 
momenta [i.e., (±27r/L, 0) and (0, ±27r/L) for the square lattice]; it is normalized to 
be comparable to its zero-momentum analogue M.^.- We define analogous observables 
for the r spins and for the composite operator ar. 

At the 4-state Potts line we have a symmetry under permutations of {a,T,aT). 
Thus, the natural choice of observables are those invariant under this symmetry. We 
have measured the following observables: 

^ I {Ml + Ml + Ml) (4.9) 

= l{S„ + Sr + S^r) (4.10) 



M' 
£ 



(4.11) 



These observables coincide with the usual ones for the 4-state Potts model up to some 
multiplicative constants. We then define the magnetic susceptibility 



the two-point correlation at the smallest nonzero momentum 

1 



V 



the second-moment correlation length 



r(2) 



1/2 



2 sin f ' 



(4.12) 



(4.13) 



(4.14) 
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the energy density (per bond) 
and the specific heat 



E = ^{£) , (4.15) 



Ch = ^ {{£') - {£)') . (4.16) 

In all these formulae, \^ = is the number of lattice sites and 2V is the number of 
bonds (we have a square lattice with periodic boundary conditions). 

Finally, let us define the quantities associated with the Monte Carlo dynamics. 
Given an observable O, we define the corresponding unnormalized autocorrelation 
function as 

Caoit) = {OsOs+t) - {Of , (4.17) 

where all the expectation values (■) are taken in equilibrium, and t is the "time" in 
units of MC steps. The associated normalized autocorrelation function is 

Poo(t) = ^ . (4.18) 

The integrated autocorrelation time for the observable O is defined as 

^int.o = IT E Poo{t) . (4.19) 

^ t=-oo 

The integrated autocorrelation time controls the statistical error in MC estimates of 
the mean (C). Given a time series of measurements {Oi, O2, ■ ■ ■ , On} (in equilibrium), 
the sample mean 

— 1 " 

O = - Y^Ot (4.20) 

^ t=i 

constitutes an unbiased estimator of (O), and its variance (when n ^ Tint.o) is 

var(O) = - 2rint,oCoo(0) . (4.21) 
n 

Thus, the variance of O is 2rint,c' larger than it would be if the measurements were 
uncorrelated. We can numerically obtain reliable estimates of both Ti^t,o and its error 
bar (from the autocorrelation function) by using a self-consistent truncation procedure 
||38| , Appendix C]. We have used a window of width QT[^t,Oy which is sufficient whenever 
the autocorrelation function decays roughly exponentially. This almost-exponential 
decay has been confirmed numerically in p4| . 

To compute the specific-heat error-bar we used the following procedure: first we 
computed the mean energy {£), and then considered the observable O = {£ — {£)Y 
using using the general procedure described in this section. 
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4.3 Summary of the simulations 



We have simulated the 4-state Potts model on an L x L square lattice with periodic 
boundary conditions, using the "embedding" algorithm described in Section 4J.. We 
performed the simulation at 194 pairs {J,L). The values of the coupling constant 
J range from 0.261 to the critical point Jc = |log3 ~ 0.274653. The lattice sizes 
L range from 16 to 1024 at the critical point, and from 32 to 512 off criticality. 
In all cases we have started our simulations with a random configuration, and we 
have discarded the first 10^ iterations to allow the system to reach equilibrium. This 
discard interval is more than sufficient: in the worst case {L = 1024 at Jc) it is roughly 
equal to 190rint,£: (or 160rexp,£; and in all other cases it is at least 300rint,£-. 

The length of the runs ranges from 9 x 10^ to 10^ iterations. For L = 512 this 
run length corresponds to 10^ times Tij^t,£] for L = 256 it is at least 2 x lO^rint,^-; for 
L = 128, at least 3 x lO'^rint,^:; and for L = 64, at least 4 x 10''rint,£-. These run lengths 
are more than sufficient to get a fairly good determination of the dynamic quantities, 
and to get high-precision data for the static quantities. Unfortunately, for L = 1024 



we were able to achieve only 1500rint,£-, as the autocorrelation time is quite large 
The data at the critical point were employed already in our previous work [p^ . 



to extract the dynamic critical behavior of our SW-type algorithm. Here we have 
improved the statistics at L = 64, 128, 256; the revised data at criticality are in 
displayed in Table |l[ The whole set of data, including the 187 runs off criticality, can 
be obtained from the authors. 

The CPU time required by our program is approximately lOL^ yus/iteration on 
an IBM RS-6000/370. The total CPU time used in this project was approximately 
8.5 years on this machine. The simulations were mainly run on the CAPC cluster 
at New York University, the IBM SP2 cluster at the Cornell Theory Center, and the 
DEC Alpha cluster at the Pittsburgh Supercomputing Center. 



5 Numerical Results 

5.1 Finite-size scaling at criticality 

In this subsection we are going to test the FSS predictions (3.9) / ( P^23D / ( ^.24]) 
for the 4-state Potts model at criticality. In p3| , |3^ we presented some preliminary 



results. However, we have now made significant extensions of these runs, in some 
cases doubling the statistics. More importantly, we now have a better theoretical 
knowledge of the FSS behavior of this model, which includes both multiplicative and 
additive logarithmic corrections. 

For each quantity O, we shall carry out fits to several different Ansatze using 
the standard weighted least-squares method. As a precaution against corrections to 
scaling, we impose a lower cutoff L > Lmin on the data points admitted in the fit, and 
we study systematically the effects of varying L^in on both the estimated parameters 
and the x^- general, our preferred fit corresponds to the smallest Lmin for which 
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the goodness of fit is reasonable (e.g., the confidence level0 is ^ 10-20%) and for 
which subsequent increases in L^m do not cause the to drop vastly more than one 
unit per degree of freedom. 

5.1.1 Second-moment correlation length 

The quantity ^^"^^ L is expected to approach a constant as L — > oo, with 
additive 0(l/logL) corrections (3.9). The constant x* can be in principle computed 
via conformal field theory, but to our knowledge this calculation has not yet been 
done. 

A fit of this ratio to a constant is reasonably good only for Lmm = 128: 

^(2) 

i_ = X* = 1.0022 =b 0.0023 (5.1) 

with = 2.38 (3 DF, level = 50%). However, if we take into account the 1/logL 
corrections we get a good fit already for Lmm = 16: 

= (1.0221 ± 0.0061) (5.2) 
L log-L 

with = 2.32 (5 DF, level = 80%). On the other hand, an equally good fit can be 
obtained using a power-law correction L~'^ with < A < 0.5; this is to be expected, 
as a logarithm can be well mimicked by a small power. We can conclude conservatively 
that 

X* = 1.02 ±0.03 . (5.3) 

5.1.2 Susceptibility 



The expected behavior of the critical susceptibility at finite L is given by (|3.24| ). 
Thus, in addition to the standard L'^^'^ term with ^jv = 7/4, we have a multiplica- 
tive logarithmic correction (logL)"^/^, and also additive logarithmic corrections of 
orders log log L/ log L and 1/logL. All these features make the accurate numerical 
estimation of the exponent 'yjv quite difficult. 

If we fit the susceptibility to the naive power law AU ^ we get a reasonably good 
fit for Lmin = 16: 

- = 1.744 ±0.001 (5.4) 

V 

with = 2.21 (5 DF, level = 82%). The difference from the exact result 1.75 is 
not large, but it is six standard deviations and thus strongly statistically significant. 
The fact that the estimate (|5.4| ) estimate is smaller than the exact value is consis- 
tent with the existence of a multiplicative logarithmic correction raised to a negative 



^ "Confidence level" is the probability that would exceed the observed value, assuming that 
the underlying statistical model is correct. An unusually low confidence level (e.g., less than 5%) 
thus suggests that the underlying statistical model is incorrect — the most likely cause of which 
would be corrections to scaling. 



22 



power. However, without the theoretical knowledge of a multiplicative logarithmic 
correction, we could equally well conclude that this small discrepancy is due to addi- 
tive corrections to scaling . 



If we try to fit the quantity x/L'^/'^ to an arbitrary power of logL [i.e., xl^^^^ = 
A(logL)^] we obtain a reasonable fit already for Lmin = 16: 

p = -0.0236 ± 0.0040 (5.5) 

with = 3.64 (5 DF, level = 60%). This result differs drastically from the theoretical 
prediction p = —1/8 = —0.125. Alternatively, we can try to fit x(logL)^/® to a power- 
law AL"'^'^; the result for Lmin = 128 is 

^ = 1.765 ±0.003 (5.6) 

with = 1.83 (2 DF, level = 40%). Once again this differs by five standard deviations 
from the exact value 1.75. 

Looking at the FSS prediction ( p.24| ), it is plausible to think that the reason why 
3D / ( ^.6| ) differ so radically from the leading-order theoretical prediction L''/^(log L)"^^ 



are the large corrections to scaling. To test this idea, we have fitted x/[L''/^(log L) ^/^] 
to A + S log log L/ log L. The fit is reasonable for Lmm = 64: 

^ (1.673 ±0.033) - (1.056 ± 0.098) ^"f^°f^ (5.7) 



L7/4(logL)-i/8 ^ • • / V • • MogL 

with = 2.76 (3 DF, level = 43%). However, the ratio B/A ^ -0.63 is quite 
different from the expected value a'/(8a^) = —1/16 = —0.0625. 

Furthermore, if we try fitting x/[L'^/^(logL)^^/^] to A + B/\ogL, we obtain an 
equally good (even slightly better) fit: with the same Lmin = 64, we get 

(1.454 ± 0.013) (5.8) 



L7/4(logL)-V8 ^ • • ^ logL 

with = 1.94 (3 DF, level = 59%). Clearly, at these modest values of L, it is 
virtually impossible to disentangle numerically the log log L/ log L and 1/logL con- 
tributions, both of which are predicted theoretically to be present. Indeed, the small 
value (—1/16) of the universal log log L/ logL coefficient makes it undetectable in the 
presence of the much larger (~ —0.4) nonuniversal 1/ logL coefficient. To distinguish 
these two contributions, we would need to reach at least log logL ^ 5, i.e. L ^ lO^''! 
Note, finally, the large discrepancy between the estimates (^.7|) and (|5.8|) of the lead- 
ing amplitude A: it is virtually impossible to estimate the correct limiting value in 
the presence of such strong corrections. 

5.1.3 Specific heat 

The FSS behavior of the specific heat at criticality is given by ( p.23| ); it is of 
the same LP(logL)~'' form as the susceptibility, but with a much larger logarithmic 
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exponent q. This makes it extremely difficult to estimate accurately the leading 
exponent a/z/ = 1, as was found already in . 

A naive power-law fit to AL'^/" gives a reasonably good result for L^m = 128: 



a 

V 



0.770 ± 0.008 



(5.9) 



with = 1.15 (2 DF, level = 57%). This estimate is far below the exact value; this 
deviation is roughly consistent with the predicted multiplicative logarithmic correc- 
tion in both sign and magnitude: a behavior (logL)^^/^ can be well mimicked over 
the interval 128 < L < 1024 by a power L-^-^^. 

We can try to estimate the power of the logarithmic term by fitting Ch/L to 
Alog^L. The result is good for Lmin = 128: 



P 



-1.258 ±0.044 



(5.10) 



with = 0.66 (2 DF, level = 72%). We obtain a much better result than for the 
susceptibility, but ( |5.10| ) is still six standard deviations away from the predicted value 
—1.5. We can also try to compute a/z/ by fitting Ch log^"^^ L to the power-law AL"^^ . 
The result for Lmm = 128 is 

(5.11) 



- = 1.044 ±0.008 

V 



with = 0.97 (2 DF, level = 62%). This is still five standard deviations away from 
the exact value a/u = 1 . 

Finally, we have fitted the function C h / [L(\og L)^^^"^] to the Ansatz A+B log log L/ logL. 
The fit is good for Lmin = 128: 



C 



H 



L(logL)-3/2 



1.607 ±0.113) - (1.926 ±0.354) 



log log L 



(5.12) 



logL 

1.2 is rather different 



with = 0.72 (2 DF, level = 70%). Again the ratio B/A s 
from the expected universal value 3a'/(2a^) = —3/4 = —0.75 from ( |3.23| ). 

If, instead, we try to fit Ciy/[L(logL)"^/^] to the Ansatz A + B/logL, the fit is 
reasonable already for Lmin = 64: 



C 



H 



L(logL)-3/2 



;i.291 ±0.022) 



1.526 ±0.096 
logL 



(5.13) 



with = 1-22 (3 DF, level = 75%). Again, we are unable to disentangle the two 
types of corrections. 

5.1.4 The combination x/cU^"^ 

From ( 3.23 ) / ( 3.24 ) we see that the combination x/C^^'^ does not show any multi- 
plicative logarithmic correction, nor any additive correction of order log log L/ log L. 
Rather, 



C 



X 

1/12 
H 



(0,0,7/;,L-1)~L5/3 



1±0| 



loffL 



(5.14) 
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It is therefore of some interest to see whether the fits for this particular combination 
are at all "cleaner" than those for x ^^nd Ch separately. To compute the error bars 
on x/C]/^^, we took into account the cross-correlations between the specific heat and 
the susceptibility. We can thus trust value in these the fits. 

To test the prediction ( |5.14| ), we first tried a power-law Ansatz AL^ . The fit is 



already good for Lmin = 16: 

p = 1.682 ±0.001 (5.15) 

with = 3.60 (5 DF, level = 61%). This exponent is 15 standard deviations away 
from the expected value 5/3 = 1.666667. This discrepancy might be due to a rather 
strong correction to scaling. 

If we assume that the corrections to scaling are of order 1/ logL, the fit is reason- 
able for Lmin = 64: 

xCrr^^^^ , ^ 0.429 ±0.057 

= (1.422 ±0.013) (5.16) 

with = 2.23 (3 DF, level = 53%). However, a correction of order log log L/ log L 
gives an even better fit: already with L^j^ = 32 we have 

yCrr^^^^ , , , , log log!/ 

^ ,f = 1.588 ±0.026 - 0.785 ±0.078 f , 5.17 
^ ' ^ logL ^ ' 

with x^ = 2.37 (4 DF, level = 67%). Again, we find it impossible to distinguish 
clearly between these two corrections to scaling, even when we know on theoretical 
grounds that only 1/logL corrections are present. 



5.2 Finite-size scaling off criticality 



In this section we are going to test the FSS predictions for the 4-state Potts model 
off criticality. In particular, from (p.ll|) we see that 



e(0,O,V^,L-i) 



3 3 
2 ' 4 



(logL; V'),0,0,1 ±0 



logL 



(5.18) 



The first argument on the right-hand side is equal (for large L) to 



0/| 3(logL;?/') 



[-a^ logL) 3/4 



a' log log L 
log L 



+ 0\ 



logL^ 



(5.19) 



This means that if we plot ^'•^V-^ a function of cpL'^^"^ log~'^''^L, the points will col- 
lapse onto a single curve. Of course, the expected deviations of order log logL/ logL 
may make this collapse less than perfect. [| 

^ One could also plot ^^^^/L as a function of 0L'^/^(logL)~'^/^[l + (a'/a^) log logL/ logL], and 
get only 1/logL corrections. However, we have seen in the preceding subsection that it is very hard 
to disentangle these two effects at criticality. Thus, we expect no gain in introducing explicitly the 
log logL/ logL correction. 
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Analogously, from ( p.21| )- (|3.22| ) we obtain similar predictions for the specific heat 
and the susceptibility: 

Ch{(P,0,iP,L-^) ~ ( \ ^/loglogL\ 

= FcA ,/ 0,0,1 +0 f ^ (5.20) 



L(logL)-3/2 

X(0,O,V',L-^) 
L7/4(logL)-V8 




(5.21) 



where and are certain derivatives of the free energy /sing- Again, plotting the 
l.h.s. of ( ^.20| ) or (|5.21| ) versus 0L^/2(logL)~^/^ will collapse the corresponding data 



onto a single curve, modulo corrections O (log log L/ log L). 

The effect of the marginal scaling field ip on the FSS equations is thus threefold: 1) 
The specific heat and the susceptibility have a multiplicative logarithmic correction. 
2) The scaled-temperature variable appearing on the r.h.s. of the FSS equations is 
not merely L}^^{J — Jc), but has a multiplicative logarithmic correction {\ogL)~^^^ 
[the same for all observables] . 3) The corrections to finite-size scaling are not 0{L~^), 
but rather 0(loglogL/ logL), which makes the numerical analysis much harder. 

The naive approach to off-criticality FSS — neglecting all multiplicative loga- 
rithms — would be to plot ^^^^/L, Ch/L, and x/L'^^'^ versus L^/'^{J — J^). These 
plots are displayed in Figure |l|. The results are quite poor: there is no data-collapse 
for the specific heat; and for the susceptibility and the correlation length the data 
collapse only close to Jc- 

A slightly less naive approach is to incorporate the predicted multiplicative loga- 
rithms for the specific heat and the susceptibility, while still ignoring the multiplicative 
logarithmic corrections to the abscissa L^^'^{J — J^,). But if we do this, we get even 
worse plots (see Figure 0), with the exception of the specific heat near Jc (where 
the rescaling of the abscissa makes no difference anyway, and the multiplicative log- 
arithm in Ch helps a lot). The slight deterioration in the susceptibility plot near 
is a reflection of our inability ( |5.5| ) to verify the correct power of log L. 

However, when we try to verify the correct FSS equations ( ^.18| )- (|5.21| ), the re- 
sult is very different (see Figure The data for the specific heat exhibit a good 
collapse away from Jc. Close to the critical point we see large deviations, presumably 
due to the log log L/ log L and 1/logL corrections. For the susceptibility, we get a 
good (though not perfect) data-collapse along the entire curve; the data-collapse is 
even better for the correlation length. We conclude that the data are in reasonable 
agreement with the predicted multiplicative logarithms, if we make allowance for the 
additive logarithmic corrections to scaling. 

Remark. In most situations, FSS plots of this type are difficult to obtain because 
of the uncertainties in the determination of Jc- a small error on Jc ruins the data- 
collapse. Fortunately, in our case we do know the exact value of Jc- This fact allows 
us to test the FSS predictions very accurately. 

Another way to present FSS data is to avoid using the "bare" variable J — Jc, and 
to use instead the physical observable ^^'^\L)/L. From ( p.29| )- p.30| ) we conclude that 
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plotting Ci^/[L(logL)~^/^] or x/[L''/^(logL)~-^/'^] versus ^^^V-^ '^ill provide a single 
curve for each observable, modulo corrections of order log log L/ log L. In Figure ^ 
we display the FSS plots of these two observables, when we neglect the multiplicative 
logarithmic corrections (left column) and with the full leading terms (right column). 
Including the multiplicative logarithmic corrections makes a big improvement, espe- 
cially for the specific heat. Again we see large corrections to scaling near the critical 
point (namely, at ■C^^V-^ ^ x* ^ 1-02), in agreement with the behavior found in 
Section pTT] . 

The main drawback of this type of plot is that data with small / L are artificially 
compressed. One way to overcome this difficulty is to replace L by ^^'^^ in the y-axis. 
That is, we can plot Cn/i^^^Klog^^^^)-^/^] and x/ [e^'^'^' (log ^^'^ versus ^^^^/L 
(see Figure H). We emphasize that these plots are identical to those in Figure ^ 
except that the scale of the y-axis is changed in a way depending on x = ^^"^^L. We 
see that the data-collapse at small ^^^V-^ is not in fact as good as it had appeared in 
Figure ^: in reality, it is mediocre for the susceptibility and quite poor for the specific 
heat. Nevertheless, there is a clear improvement in both cases when we include the 
predicted multiplicative logarithmic correction. 



5.3 Extrapolation techniques 

We have also considered how well the extrapolation scheme introduced in jlD 



works in the presence of logarithmic corrections. This method makes use of the 
following FSS equation for an arbitrary long-distance observable O: 



0{J,2L) ^ (i{J,L)- 



0{J,L) ^""y L 



+ 0(r",^""), (5.22) 



where Fo is an unknown scaling function and is a correction-to-scaling exponent. 
(Here we have chosen a size-scaling factor s = 2; a similar equation holds, of course, 
for any s.) 

Let us first check the analogous equations for our case. To simplify the notation, 
we will write only two arguments for the observables: 0{J, L) = 0{(j), h = 0,iIj, L~^). 
From ( 3.17 ) we conclude that 

where and are universal (albeit unknown) scaling functions. Using the relation 
( p.l8| ) we arrive at 

where and are again universal. Thus, we get the same equation as ( [5. 221 ), but 
the corrections to scaling are much larger (~ 1/logL). 

The same type of equation can be derived for the specific heat and the suscep- 
tibility, using (|3.29| )- (|3.3CI| ): although there are corrections of order loglogL/logL, 
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these cancel out in forming the ratio between lattice sizes L and 2L. Notice that the 
nonuniversal /co/logL contribution also cancels out in the ratio 0{J,2L)/0{J, L). 
This means that the leading correction to scaling is universal. That is, we have 

__ = .4o^— j+Bo^— j — + ... (5.25) 

where both scaling functions Aq and Bq are universal. 

The 1 / log L corrections in ( p.25| ) make a numerical determination of the finite- 
size-scaling functions Ao very difficult. Unfortunately, these functions are the basic 
objects needed to extrapolate from finite L to the infinite- volume limit [|l9l. Thus, a 



naive application of the method of Ref. |T9| is not expected to work well. Figure ^ 
shows the plots of 0{J,2L) /0{J, L) versus ^*-^V-^ ^'^^ ^^e specific heat, the suscepti- 
bility and the second-moment correlation length. The data collapse is notably inferior 



to that observed in most other models [19|. Note also that the specific-heat curve 



is very far from its correct value 2°'^'^ = 2 at the critical point ^'•^V-^ = x* ~ 1.02; 
rather, it takes the value 2'^"/'^^<=*f = 2^'^'™ ^ 1.71. This is another indication of the 
very strong corrections to scaling. 



6 Dynamic critical behavior of the Swendsen- Wang- 
type algorithm 

In this section we shall analyze briefiy the dynamic critical behavior of our Swendsen- 
Wang-type algorithm for the two-dimensional 4-state Potts model. Recall that we 
have used the "embedding" variant of the SW-type algorithm for the Ashkin-Teller 



model ||33|, This algorithm does not coincide with the standard Swendsen- Wang 



37| algorithm for the 4-state Potts model, but it is expected heuristically (and con- 



firmed numerically ||3^ at criticality) to lie in the same dynamic universality class. 
We therefore expect the dynamic critical exponents to be identical for the two algo- 
rithms; and we also expect the dynamic finite-size-scaling functions to be identical 
modulo a multiplicative factor. 

6.1 Dynamic finite-size scaling at criticality 

We first fit the integrated autocorrelation time for the energy, rint,^:, to a pure 
power law Tint,£ = AL^'"*-^. We obtain a good fit for Lmm = 32: 

Zint,£ = 0.876 ±0.011 (6.1) 

with = 2.54 (4 DF, level = 64%). However, this cannot be the true asymptotic 



behavior, because the Li-Sokal bound |3^, |3^ guarantees that Ti^it,s ^ const x Ch, and 



we know from ( ^.23] ) that the specific heat diverges at criticality like L(logL) '^/^.Q 



^ The Li-Sokal bound has been proven for the original Swendsen- Wang algorithm ]39|] and more 
generally for the direct form of the Ashkin-Teller SW algorithm |Q. But the embedding form of 
the AT algorithm is believed (and observed numerically ||3^ ) to be in the same dynamic universality 
class as the direct algorithm. So the bound should apply for it as well. 
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Rather, this apparent exponent Zint,£ < 1 is probably an effect of multiphcative loga- 
rithmic corrections (with a negative exponent), just as it is for the specific heat. 
In p3, HOl we showed that there are only two likely behaviors for the autocorrela- 



tion time of the Swendsen-Wang algorithm for the 2D Potts models: 

, „ f ALP with p small ,„ „x 

With the available statistics it is very hard to distinguish between these two scenarios. 
If we fit the quantity Tint,£/[L(}ogL)~^/'^] to a pure power law AL^, we obtain a 
reasonable fit for Lmin = 128: 

p = 0.153 ±0.028 (6.3) 

with = 1-30 (2 DF, level = 52%). On the other hand, if we fit the same quantity 
to the logarithmic Ansatz A + BlogL, the fit is very good already for Lmin = 16: 

T-int,f (0.197 ± 0.173) + (1.278 ±0.044) logL (6.4) 



L(logL)-3/2 



with = 1-53 (5 DF, level = 77%). Thus, our MC data supports slightly better the 
logarithmic scenario compared to the power-law scenario. 

One can also study directly the ratio T[^t,£/CH, and fit the results to the Ansatze 
(3.2). Unfortunately, we do not know the correct error bar on this ratio, as we do not 



know the covariance between the estimator of the specific heat and the estimator of 
rint,^:. Instead, we shall use the upper bound on the error bar provided by the triangle 
inequality. Our error bars are thus overestimated — by how much we do not know 
— and, consequently, the values of the fits are expected to be artificially small. 
We can, however, compare the relative values for the two Ansatze. (See for 
a complete discussion concerning these points.) If we try to fit T^^t^elCu to a pure 
power law AL^ ^ we get a good result for Lmin = 16: 

p = 0.119 ±0.011 (6.5) 

with x^ = 1.06 (5 DF, level = 96%). With the logarithmic Ansatz A ± SlogL, the 
fit is less good: for the same Lmin we obtain = 1.99 (5 DF, level = 85%). Thus, 
for Lmin = 16 the power-law fit gives a a factor of two better than the logarithmic 
fit. However, for Lmin = 32, the two fits give more nearly equal values of x"^'- 0-98 
and 1.27, respectively. 



6.2 Dynamic finite-size scaling off criticality 

The FSS behavior of the dynamic observables off criticality is expected to be 
similar to that discussed previously for the static observables. That is, we expect 
multiplicative logarithmic corrections, as well as additive logarithmic corrections to 
scaling. In this section we have considered the two scenarios discussed in the previous 
subsection, with the aim of distinguishing between them. 
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First, we have made the FSS plots analogous to those of Figures and 0. That 
is, we have plotted the leading behavior — Ti^t,s/ L^"^^'^ or rint,£:/[-Z^(logL)~^/^(l + 
6.501 logL)] — versus L^^'^{J — Jq) (left column in Figure |^. In both cases, we 
obtain terrible fits. However, the situation improves dramatically when we consider 
the correct abscissa L^/^(logL)~^/^(J— Jc") (right column in Figure]^). Unfortunately, 
from these plots corresponding to the two most likely scenarios, it is impossible to 
tell which one is more likely. 

Alternatively, we can plot the same quantities as functions of the physical observ- 
able ^*^^^(L)/L (left column in Figure]^). Again, the plots are very similar, and there 
is no objective reason to prefer one over the other. The corrections to scaling for 
small ^^^V-^ appear to be quite small; they grow slightly in size as we approach larger 
values of ^^"^^ /L. To check that this behavior is not an artifact of the compression 
of the vertical axis at small ^^"^^/L, we have tried replacing L by ^'■^^ in the y-axis 
denominators (right column in Figure]^). We emphasize that these plots are identical 
to those in the left column, except that the scale of the y-axis is changed in a way 
depending on x = ^^"^y L. Contrary to what we observed for the static observables, 
the corrections to scaling for small C*^^V-^ large compared to the error bars, 

except for a few points at very small in the logarithmic Ansatz (where the neglect 
of the additive constant A in A + Blog^^'^^ obviously plays a role). Note, finally, that 
the scaling function tends to a nonzero constant when ^'•^V-^ ~^ 0; expected on 
general theoretical grounds. 



A Calculation of the coefficient a' 



In this appendix we derive the value ( p.lO|) of the cubic coefficient a' in the RG 



flow for the 4-state Potts model. The method is identical to that of Cardy et al. 
0, carried to one higher order: namely, we study the renormalization-group flow for 
the g-state dilute Potts model in a neighborhood of the multicritical point g = 4, 
(j) = h = ijj = 0] we then match this flow with the exactly known values of the 
critical exponents as a function of g, expanded now through second order in (4 — g)^/^. 



A.l Exactly known values of critical exponents 

The leading thermal exponent yT,i, the next-to-leading thermal exponent yT,2, 
the leading magnetic exponent yH,i and the next-to-leading magnetic exponent yH,2 
of the g-state dilute Potts model are known exactly, both for the ordinary critical 
point (arising in the pure Potts model) and for the tricritical point (arising in the 
dilute model). Their values are: 

yr, = ^ (A.1) 

4a; , . , 

VT, = (A.2) 
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15 + 8x + x2 

^-■^ = 4(2 + x) (^-'^ 

^^^.2 = , (A.4) 

where 

x = — arccosj^^) , (A. 5) 

with — 1 < a; < corresponding to the ordinary critical points and < x < 1 corre- 
sponding to the tricritical points. The value of yT,i for the ordinary critical point was 
conjectured by den Nijs |0; this conjecture was extended to the tricritical branch of 
the dilute Potts model by Nienhuis et al. 0] . It has now been derived by Coulomb-gas 
methods pil, H2| for both the critical P3 and tricritical E3| branches. The expression 




( |A.2| ) for the next-to-leading thermal exponent was first conjectured by Burkhardt 
||43|| ; it was derived using Coulomb-gas methods by Nienhuis [^]. Finally, the ex- 
pressions ([A.3|) / ( [A~4D for the leading and next-to-leading magnetic exponents were 
first conjectured by Nienhuis et al. and Pearson [||; their validity was established 



by den Nijs [0] using the Coulomb-gas approach. In the Coulomb-gas formulation, 
all these formulae are naturally parametrized in terms of the variable t = 1/(2 + x), 
which is proportional to the temperature of the associated Gaussian model; the ther- 
mal (resp. magnetic) exponents are linear (resp. quadratic) polynomials in t. 

These exponents can also be understood in terms of conformal field theory |^ 



46| , 0. Let us consider a conformal field theory (CFT) with central charge c < 1, 
which we parametrize as 

c = 1 (A.6) 

m[m + 1) 

with < m < cxD (not necessarily integer). Then the conformal weight of a primary 
field {t, s integer) is given by the Kac formula^] 

_ [{m + l)r - msf ~ I 
Am[m + Ij 

The corresponding critical exponent is y^.s = 2 — 2A,f. ,,. If m is a rational (resp. 
irrational) number, then the corresponding CFT has a finite (resp. infinite) operator 
algebra. If m is an integer > 2, then the corresponding QFT is unitary, at least as 
long as we restrict attention to the fields 0^,8 satisfying l<s<r<m — 1. 

The critical g-state Potts model can be represented |^8|, ^ by a CFT with central 
charge 

c^l-^. (A.8) 



where the parameter x is given by ( |A.5|) with — 1 < x < 0.0 We make the identification 



X = -2/(m + 1) , (A.9) 



^ We follow the notation of |46| ^ . In Ref. |^ the indices r and s are interchanged with respect 
to this notation. 
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Our parameter x is equal to minus the parameter y appearing in 
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and the Kac formula ( |A.7| ) reduces to 

_ [2{T-s) + sxf-x^ 
^'■'^ - 8(2 + x) ■ ^^-^^^ 

The unitary theories m = 2, 3, 4, 5, . . . correspond to g = 1, 2, 4 cos^(7r/5), 3, . . . . 

The energy operator e is identified with 02,i- This means that yx.i = 2 — 2A2,i, 
and we obtain the resuh ([A.l| ). From the conformal algebra ee ~ 02,i02,i ~ / + 



03^1, and we identify the operator (p^^i as the one giving the next-to- leading thermal 
corrections. This implies that yT,2 = 2 — 2A31, and we get ( [A .21 ). Finally, the 
spin operator s is identified with 0m+i m+i (or alternatively ). Therefore, 

Uh,! = 2 — 2A m+i m+i = 2 — 2A m+i m+i_-^, aud this yields ( |A.3D . The second magnetic 
exponent yH,2 can be derived from the conformal weight of the operator 
(or alternatively ). More generally, we can obtain the full spectrum of 

thermal exponents 

(4 - n^) + (n + 2)x 
yT,n = 2-2A„+i,i = ^ 1—^ (A.ll) 

and magnetic exponents 

[16- (2n - 1)2] + 8x + x2 



|///,n = 2 - 2A™ii+„_i_m±i = 2 - 2Am±i_„_m±i 4(2 + a;) 

(A.12) 

However, this identification of the magnetic operators makes sense only when m is 
an odd integer 1, 3, 5, 7, ... , corresponding to g = 0, 2, 3, 2 + y^, .... 

The tricritical Potts model can be obtained by the analytic continuation of the 
preceding formulae to < x < 1 The conformal charge continues to be given by 
but we must now identify 



X = 2/m . (A.13) 
The Kac formula for the tricritical models is thus given by 

_ [2{r ~ s) + rxf - x^ 

^^'^ - W+^) ■ ^ ^ ^ 

The energy operator is now identified with the operator 0i 2, which reproduces (|A.1|) . 
More generally, the thermal exponents are given by the operators 0i^i+„, while the 
magnetic exponents come from the operators 0™^™+„ (or equivalently from 0™ _™+i_„). 
Notice that the magnetic operators make sense only when m is an even integer 
2, 4, 6, 8, ... , corresponding to the same series q = 0,2,3,2 + ^2,.. . as in the critical 
case. 



A. 2 Renormalizat ion-group flow for the g-state dilute Potts 
model 

Following Cardy et al. H, Section II], we study the RG fiow for the g-state dilute 
Potts model in a neighborhood of the multicritical point e = g— 4 = 0, = /i = ^/' = O. 
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Cardy et al. showed that the RG equations can be brought by a series of smooth 
changes of variable into the form 

dip 



dl 



a{ilj^ + e) + 0(f ,e/,e2) (A.15a) 
{yT + lnP)<P + 0(f,e/,e2) (A.15b) 
{yH + cilj)h + 0(/3,e/,e2) (A.15c) 



dl 
dh 

'dl 

where a, b, c, yx, yu are universal parameters. Here / stands collectively for the fields 
and h] and we consider e to be of order ip'^ because that is its magnitude at the 

fixed points = — e + O(e^). 

By a similar argument one can treat the terms of cubic order in the fields; as 

before, the if) equation decouples from the others, and becomes 

^ = a(V>2 + e) + + a'e^j + 0(V^e^/.^e2) . (A.16) 
For e < there are fixed points at 



^± = H-e?" + ^(-^) + Oii-ef/') (A.17) 
with eigenvalues 

y± = ±2a(-e)i/2 + 2a'(-e) + 0{{-ef/^) . (A.18) 

Here the + sign corresponds to the tricritical point, with a relevant second thermal 
exponent y+ > controlling the tricritical-to-critical crossover; and the — sign cor- 
responds to the ordinary critical point, with an irrelevant second thermal exponent 
y_ < controlling the leading corrections to scaling. Matching ( [A.18|) to the known 
second thermal exponent ( |A.2| ) / ( |A.5| ) , we find 

a = - (A. 19) 

TT 

a' = --^ (A.20) 

This reasoning gives an independent (and consistent!) derivation of the quadratic 
coefficient a, which was derived previously by Nauenberg and Scalapino 0, ^ by 
matching ( [A.15| a,b) to Baxter's ^ exact result for the latent heat in the pure Potts 
model for g > 4; and it gives the promised derivation of a'. 

Remark. This reasoning does not fix the value of a". Indeed, by a smooth 
change of variable ip = ip' + «2V''^) the coefficient a" can be set to any desired value: 
for example, a" = or a" = a'. The latter choice has the property of placing the 
fixed points at ip± = ±(— e)^/^ with no correction 0(— e). This approach can also be 
extended to higher order: e.g. starting from the fourth-order equation 

^ = ai^p"^ + e) + a'tp^ + a"eip 

+ o(3)^4 ^ ^(4) ^^2 ^ ^(5)^2 ^ ^^^5^ ^^3^ ^2^^ ^3) ^^^21) 
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we can make a smooth change of variable ip = {1 + aue)ip' + 02'^'^ + cts"^'^ and fix 
three of the four coefficients a" , a'^^\ a^^\ a^^^ to any desired values. For example, if 
we choose to fix a" = a', a^^^ = a^'^\ a^^^ = we then have the normal form 

^ = + e)F(^) + 0(^^ eV, e') (A.22) 

with F independent of e. Alternatively, we can choose a" = a', a^^^ = 0, a^^^ = a*-^-* to 
yield the normal form 

^ = (^2^e)[Ao(6) + Ai(e)^] + 0(^^e^^e^^,e=') (A.23) 

with y4o, Ai independent of ip. It would be interesting to know whether either (or 
both) of these normal forms can be established to all orders in and e. If so, the 
functions F{iIj), Ao{e) and Ai{e) would then be completely determined from the known 
exponent (^/(^. 
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L 


MCS 


A 


Cii 




' iiii 


16 


0.9 


141.41 ± 


0.29 


2.513 ±0.013 


15.758 ± 


0.056 


12.86 ± 


0.24 


32 


1.9 


474.23 ± 


0.94 


4.170 ±0.020 


31.681 ± 


0.100 


23.13 ± 


0.40 


64 


9.9 


1587.70 ± 


1.90 


6.971 ±0.020 


63.704 ± 


0.115 


41.58 ± 


0.42 


128 


6.9 


5321.74 ± 


10.55 


11.730 ±0.055 


128.153 ± 


0.368 


77.13 ± 


1.26 


256 


4.4 


17764.70 ± 


60.96 


19.936 ±0.158 


256.291 ± 


1.243 


143.00 ± 


3.99 


512 


2.9 


59876.54 ± 


334.67 


33.969 ±0.467 


519.102 ± 


4.078 


252.83 ±11.57 


1021 


0.8 


196872.13 ±3137.60 


60.135 ±2.091 


1020.767 ±21.680 


531.11 ±67.68 



Table 1: Monte Carlo data for the 4-state Potts model at criticality. For each lattice 
size L we show the number of measurements (MCS) in units of 10^, the susceptibihty 
X, the specific heat Ch, the second-moment correlation length and the integrated 
autocorrelation time for the energy Tint,f . The quoted error bars correspond to one 
standard deviation (i.e., confidence level ~ 68%). 
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d=2, q=4 Potts model: Specific Heat 
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d=2, q=4 Potts model: Correlation Length 



1.0 



> 0.5 H 



0.0 



01 ^ + 

® a, + 



-3 



-2 -1 
L^/^ (J-JJ 



Figure 1: Naive finite-size scaling plots, neglecting all multiplicative logarithms, for 

the specific heat Ch, the susceptibility x? and the second-moment correlation length 
Note that he abscissa L^/'^{J — Jc), where Jc= \ log 3 = 0.274653 . . . , also lacks 
the predicted multiplicative logarithms. Symbols denote the different lattice sizes: 
L = 32 (+), L = 64 (x), L = 128 (□), L = 256 (O), and L = 512 (o). 
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d=3, q=4 Potts model: Specific lieat 
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Figure 2: Finite-size scaling plots that include the predicted multiplicative loga- 
rithms Ch ~ L(logL)~^/^ and x ~ L'^^^(\ogL)~^^^ but do not include any multiplica- 
tive logarithms in the abscissa L^^^{J — Jc). Symbols denote the different lattice sizes: 
L = 32 (+), L = 64 (x), L = 128 (□), L = 256 (O), and L = 512 (o). Note that the 
fits are worse than in Figure 0. 
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d=2, q=4 Potts model: Specific Heat 
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d=2, q=4 Potts model: Susceptibility 
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d=3, q=4 Potts model: Correlation Length 
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Figure 3: Finite-size scaling plots for the specific heat Ch, the susceptibihty x, and 
the second-moment correlation length ■^'^^^ that take into account all predicted multi- 
plicative logarithms. Note that the correct abscissa variable is {J — Jc) L^^'^ {log L)~^/^ . 
Symbols denote the different lattice sizes: L = 32 L = 64 (x), L = 128 (□), 
L = 256 (O), and L = 512 (o). 
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d=S, q=4 Potts model: Specific Heat 
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Figure 4: Finite-size scaling plots for the specific heat Ch and the susceptibility x 
as functions of ^'^'^\L)/L. For each observable, the plot in the left column includes 
only the leading power law, while the plot in the right column takes into account the 
predicted multiplicative logarithmic correction. The improvement is clear. Symbols 
denote the different lattice sizes: L = 32 (+), L = 64 (x), L = 128 (□), L = 256 
(O), and L = 512 (o). 
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d=2, q=4 Potts model: Specific Heat 
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d = 2, q=4 Potts model: Susceptibility 
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Figure 5: Finite-size scaling plots for the specific heat Ch and the susceptibility x as 
functions of ^^'^\L)/L. These are the same plots as in Figure ^ with the lattice size 
L replaced by the second-moment correlation length ^(^^(L) on the y-axis. Symbols 
denote the different lattice sizes: L = 32 (+), L = 64 (x), L = 128 (□), L = 256 
(O), and L = 512 (o). 
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d=2, q=4 Potts model: Specific Heat d=2, q=4 Potts model: Susceptibility 
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Figure 6: "L/2L" finite-size scaling plots, of the kind employed in Ref. [0, for the 
specific heat Ch, the susceptibility x ^^^^ the second-moment correlation length 
For each observable O, we plot the quantity 0{2L) /0{L) versus ^^^^L. Symbols 
denote the different lattice sizes: L = 32 (+), L = 64 (x), L = 128 (□), L = 256 
(O), and L = 512 (o). 
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Autocorrelation time t=_j. t-: z = 0.876 
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Autocorrelation time Tj^^^j.: Log Ansatz 
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Figure 7: Finite-size scaling plots for the integrated autocorrelation time Tint,£ as 
a function of L^/'^{J — Jc) (left column) and (J — Jc)L^^^(\.ogL)~^/'^ (right column). 
We show the two most hkely scenarios: Tint,£ ~ L^'"*>^ with Zint,£ ~ 0.876 ± 0.011 
(upper row), and Tint,s ~ -L(logL)~^/^(l + 6.501 logL) (lower row). Symbols denote 
the different lattice sizes: L = 32 (+), L = 64 (x), L = 128 (□), L = 256 (O), and 
L = 512 (o). 
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Figure 8: Finite-size scaling plots for the integrated autocorrelation time Ti-at,£ as 
a function of the physical observable ^^'^\L)/L. The upper row shows the Ansatz 
Tmi,,£ ~ L^'"*'^ with -Zint,^ = 0.876 ± 0.011, while the lower row shows the scenario 
Tint,£' ~ L(logL)~^/^(l + 6.501 logL). In the left column we show the usual FSS plots; 
in the right column we have substituted the variable L in the y-axis denominator by 
in order to amplify the region with smaller ^"-^V-^- Symbols denote the different 
lattice sizes: L = 32 (+), L = 64 (x), L = 128 (□), L = 256 (O), and L = 512 (o). 
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